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^ ■ Using the measured spectrum shape for B — > tt£z/, the rate for 5+ — > 

^ . tt+tt , information on the Cabibbo-Kobayashi-Maskawa (CKM) matrix ele- 

j ment |V^|, and theoretical inputs from factorization and lattice gauge theory, 

we obtain an improved estimate of the "tree" contribution to B° — > tt+tt - . 
t-h I We find the branching ratio B(B° -> tt+tt-) | tree = (5.25^) x 10 -6 , to be 

^ ■ compared with the experimental value B(B° — ► tt+tt - ) = (4.55±0.44) x 10 -6 . 

■ The fit implies = (3.62 ± 0.34) x 10 -3 . Implications for tree-penguin 

interference in B° — > tt + tt~ and for other charmless B decays are discussed. 

in 
o . 
m ■ 

PACS Categories: 13.25.Hw, 14.40.Nd, 14.65.Fy, 11.30.Er 
^ I. INTRODUCTION 

Oh: 

The semileptonic process i? — > tt£z/ involves a form factor F + (q 2 ) related for q 2 = -m\ 
to the factorized color-favored "tree" contribution in 5° — > tt + 7t _ 13 E] • in previous 
work jl] we obtained an estimate of this contribution implying a branching ratio B(B° — > 
7r+7r ~)|tree = (7.3 ± 3.2) x 10~ 6 . A measurement of the spectrum dT(B — > Tc£u)/dq 2 has 
now been presented by the CLEO Collaboration j3] working at the Cornell Electron 
Storage Ring. Further results are expected from the BaBar and Belle Collaborations at 
asymmetric e + e~ colliders. Using the CLEO measurement and other inputs, we find in 
the present paper an improved estimate implying B(B° — ► Tr + Tr-)| tr cc = (^>-^-o'Io) x 
10~ 6 , to be compared with the observed branching ratio B(B° — > Tr + Tr _ ) = (4.55 ± 
0.44) x 10 -6 . This result has a number of implications for tree-penguin interference in 
B° — ► tt+tt - and for other charmless B decays. 

We review theoretical inputs, including constraints from factorization and lattice 
gauge theory calculations, in Sec. II, while data are discussed in Sec. III. We perform 
a global fit to these inputs in Sec. IV. The consequences of this fit are discussed for 
B° — > tt+tt - and other charmless B decays in Sec. V. We conclude in Sec. VI. 
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II. THEORETICAL INPUTS 



The B — ► 7r matrix element is parametrized by two independent form factors: 
(7r(p)|« 7M 6|S(p + ?)> = fap + q-q™*^) F + (g 2 ) + g,^i^Fo(g 2 ) , 




For massless leptons (assumed here), only F+(q 2 ) contributes to the differential decay 
rate 

AV 11/ I 2 

^(B- - = ^L| AP |F +tf )|« , (2) 

where is the relevant CKM matrix element. In the factorization hypothesis, one 
replaces the lepton pair with a pion, giving what we term the "tree" contribution T (in 
the notation of 6j) to the nonleptonic decay B° — > 7r + 7r - . In the limit of small m n , the 
two processes are related by j2] 



r tree (5° -> 7T+7T-) = 67r 2 /'|K d r|ai 



a , 2lT , |2| _ |2 dT(B° - tt-^+i/i) 



dg 2 



(3) 



where |ai| is a QCD correction which we shall take equal to 1. The majority of QCD 
effects are expected to be associated with the form factors F + ^(q 2 ) and thus are taken 
into account by the factorization ansatz. 

Other contributions to charmless strangeness-preserving B decays which we shall 
consider include color-suppressed tree (C) and penguin (P) amplitudes. The corre- 
sponding strangness-changing amplitudes are denoted by primes. We neglect smaller 
amplitudes which involve spectator quarks. For these and other details, see, e.g., Ref. 
jSj. The amplitude for B° — > 7r + 7r~ is then 

A(B° ^tt+tt-) = -(T + P) = -|T|e 17 - \P\e- ip e iS , (4) 

where we have introduced phases of CKM elements, assuming the phase of the b — > d 
penguin to be dominated by the top quark, and 5 denotes a relative strong phase. A 
question which has been of interest for some time |H |H] is whether the small branching 
ratio for B° — > 7r + 7r~ reflects the effect of destructive tree-penguin interference. If so, 
by combining this information with CP- violating asymmetries in B° — > 7r + 7r~, one can 
learn a good deal about both weak (i.e., CKM) and strong phases 0. 

We use notation in which the square of an amplitude directly gives a B° branching 
ratio in units of 10^ 6 . The observed branching ratio B(B° — > 7r + 7r~) = (4.55±0.44)xl0 -6 
(see Sec. Ill) then corresponds to \T + P\ = 2.13 ± 0.10 in our units. In previous work 
0] we found \T\ = 2.7 ± 0.6, too large an error to display any possible tree-penguin 
interference. 

Another amplitude which will be of use to us is A(B + — > 7r + 7r°) = — (T + C) / v2. The 
color-suppressed amplitude C is expected to have small phase and magnitude relative 
to T [TU]. We shall use only the conservative range [TT] 0.08 < \C/T\ < 0.37. This will 
provide a useful bound on T based on B(B + — > 7r + 7T°). 
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Lattice gauge theories predict not only the shape, but also the normalization, of 
the B — > 7r form factors at large q 2 or small pion recoil momentum in the B rest 
frame [T2"l IT3~1 ITU ITS] . These predictions turn out to be very helpful in constraining 
parameters on the basis of the q 2 spectrum in B — > 7r/z/. However, they do not address 
the key question of the form factor behavior at small q 2 or large pion recoil. 

The CKM parameter V u b is another key input whose determination is for the moment 
still subject to theoretical uncertainties. Good understanding of the B — > n form factor 
would reduce these uncertainties. Independently of B — > ttIv (or the more complex 
process B — > plu), however, various inclusive methods have been employed to extract V u b 
from semileptonic b — > u decays, including the study of leptons with energy exceeding the 
endpoint for b — > civ, the rejection of events with recoil mass above charm threshold, and 
the use of the photon energy distribution in b — > 57 to measure the "Fermi distribution" 
of b quarks inside a B meson. These are summarized in a subsection of the Review of 
Particle Physics [To] . 

The form factor F + (q 2 ) is expected to have a B* pole, as well as possible higher-lying 
poles in q 2 . In Ref. jl] we approximated it with a dipole form proposed by Becirevic and 
Kaidalov [T7j on the basis of lattice gauge theory calculations: 



F (a 2 ) = c B (l-a B ) 
+ {q ] (l-q 2 /m 2 B ,)(l-a B q 2 /m 2 B ,) ' l ° J 

A value of a B between and 1 would correspond to a pole lying above m B *. However, 
we were unable to achieve a good fit to the CLEO B — > ttIv spectrum with this form. 
(The x 2 °f the fit is more than 3 for one degree of freedom.) A generalization of the 
above form factor is to multiply it by (1 + aq 2 /m 2 B ,), where a is an additional parameter. 
The resulting form factor is equivalent to an explicit dipole: 

l — q z /m B , ± — a B q z /m B * 

However, we were unable to achieve any fit for any physical ag to the numerical inputs 
in Sec. Ill which represented any improvement over the single pole with this form. We 
thus choose instead the two-parameter form 

F + {q 2 )= : F f) 2 (1 W/™|.) • (7) 
1 — q 1 m B « 



III. NUMERICAL INPUTS 

We summarize some information on B — > tttt branching ratios fHl E3 ED] in Table 
I. The central value of the 7r + 7r° branching ratio exceeds that of 7r + 7r~ despite the fact 
that the coefficient of its dominant T term is divided by v2 (see Sec. II). To extract an 
amplitude for comparison with B° decays we must first divide all B + branching ratios 
by the B + /B° lifetime ratio EH r+/r = 1.073 ± 0.014 . 

Our estimate of \T\ based on B + — > 7r + 7r° then proceeds as follows. After correcting 
for the lifetime ratio, we find \T + C\ = 3. 13 ±0.24. With [TT] \T + C\ = \T\ (1.23 ±0.15) 
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Table I: Branching ratios for some charmless two-body B decays, in units of 10 6 . 



Mode 


BaBar [18 


Belle [19 


CLEO [20J 


Average 


7T + 7T~ 
7T+7T 


4.7 ±0.6 ±0.2 
5.5 ± 1.0 ±0.6 


4.4 ±0.6 ±0.3 
5.3 ±1.3 ±0.5 


a C+I.4+U.5 
2-0.4 
a fi +1.8+0.6 
4 -0-l. 6-0.7 


4.55 ± 0.44 
5.27 ±0.79 



we then obtain |T| = 2.55 ±0.37. This is consistent with our previous determination [I] 
but with smaller errors. (The estimate \T\ = 3.0 ± 0.3 of Ref. [22] uses too restrictive a 
value of \C/T\ in our opinion.) We seek further information from the B — > ttIu spectrum 
shape and other sources. This value would imply B(B° — > 7r + 7r~)| t rcc = (6.5±1.9) x 10~ 6 , 
only la above the experimental branching ratio. 

The CLEO Collaboration jHj measured dB(B — > -nlv)jdq 2 in three q 2 bins, each about 
8 GeV/c 2 wide. The results are not very sensitive to the choice of form factor and we 
quote them for the form factor [23] which appears to fit the data best: 

m ( (0.431 ± 0.106) x 10~ 4 (0 < q 2 < 8 GeV 2 ) 

f dq 2 — (B° -> 7T-/ + ^) = I (0.651 ± 0.105) x 10~ 4 (8 < q 2 < 16 GeV 2 ) (8) 
J dq { (0.245 ± 0.094) x 10" 4 (16 GeV 2 < q 2 ) 

Only statistical errors (dominant) are shown. These sum to a branching ratio of B(B° — > 
w-l + ui) = (1.33 ± 0.18 ± 0.11 ± 0.01 ± 0.07) x 10~ 4 , where the errors are statistical, 
experimental systematic, pion form factor uncertainty, and p form factor uncertainty. 

Lattice calculations of the form factor F + (q 2 ) have been presented in the past few 
years by the UKQCD [TJ, APE [T3], Fermilab Q3] and JLQCD |15| Collaborations. 
Numerical values of F + (q 2 ) are computed in the range 13.6 GeV 2 < q 2 < 23.4 GeV 2 . 
Although small variations are present among the four different calculations, all results 
are consistent with each other within errors. We will include them all in our fits. 

For \V u b\ we combine determinations presented in Ref. jTHj in the following manner. 
All numbers will be quoted in units of 10~ 3 . The inclusive LEP average is 4.09 ± 0.37 ± 
0.44 ± 0.34 while the inclusive CLEO value is 4.12 ± 0.34 ± 0.44 ± 0.33, where the errors 
are statistical, experimental systematic, b —>■ c uncertainty, and b —>■ u uncertainty. In 
addition there are theoretical uncertainties estimated to range up to 15%. Combining 
the two inclusive numbers before folding in the theoretical uncertainties, and treating 
the last two errors as common, we obtain 4.11 ± 0.61. We shall use this in our fits. We 
do not include some preliminary results presented by CLEO |21] and Belle |25j . 

An earlier CLEO exclusive determination of \V u b\ utilizes both tcIu and plv decays 
[TSj . Its result, which we do not use in the present fit, amounts to an average of 3.25 with 
experimental and theoretical errors comparable to those in the inclusive determinations. 
Averaging it with the inclusive value noted above, we should expect a global fit to give 
10 3 \V u b\ — 3.68 ±0.43 with an additional 15% theoretical error, or approximately a 20% 
error overall. We shall see that a modest improvement upon this error is possible, while 
the central value does not change much. 
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Table II: Sources of \ 2 m global fit to F + (q 2 ). 



Source 


x 2 


Reference 


B — > nil/ spectrum 


2.42 





Inclusive \V u b\ values 


0.65 


PS 


UKQCD lattice data 


0.74 


C3 


APE lattice data 


0.12 


PS 


Fermilab lattice data 


4.53 


M 


JLQCD lattice data 


0.21 


.15 



IV. GLOBAL FIT 

We perform an overall three-parameter \ 2 fit to the above-mentioned B — > 7r/z/ 
branching fractions in the three bins, the averaged inclusive \V u b\, and 26 lattice data 
points on F + (q 2 ). We neglect the small correlations among the three branching fractions. 
The quality of the fit is fairly good, with x 2 = 8.7 for 27 degrees of freedom. The x 2 ' s 
contributed by specific sources are summarized in Table |H] More than 50% of the \ 2 
comes from the Fermilab lattice points, which appear to be of a somewhat different 
pattern from the other three lattice determinations. 

The results of the fit are 

1-14182 , (9) 
0.23 ±0.04 , (10) 

(3.62 ±0.34) x 10~ 3 , (11) 

(5.25 ± 1.67) x 10" 6 , (12) 

2.29 ±0.36 . (13) 

A theoretical error of ~ 15% must be added to |T4&|- The value of |T| overlaps that 
(|T| = 2.55 ± 0.37) obtained in Sec. Ill from B + — > 7r + 7r°, but the 7r + 7r° value indicates 
that B(B — > 7r + 7r _ )|trcc is no smaller than 4.75 x 10~ 6 . Hence we shall truncate our 
parameter space at this lower limit, and quote 

B(£°^7r+7r-)| tree = (5.251^) x 10- 6 , (14) 
\T\ = 2.29t°oit ■ (15) 

The ranges of parameters contributing to the global fit are illustrated in Fig. ^ where we 
show the points corresponding to minimum \ 2 = 8-7 an d the ellipses corresponding to 
Ax 2 = 1. The various projections are helpful in visualizing the full range of parameter 
variation. In particular, the value of B(B° — > 7r + 7r~) | trce can vary substantially as a result 
of the uncertainty in F(0), which is still not well constrained by the data. However, its 
la upper limit of 6.92 x 10~ 6 is well below that implied by the previous estimate of Ref. 

Ill 



F(0) = 

\V u b\ = 
B(B° ^ n+n-)\ tme = 
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Figure 1: Left column: projections of error ellipsoid for global fit on the plane of two 
parameters for central values of the third. Right column: ellipses involving the tree 
amplitude corresponding to the variations shown in the left column. Note that |T| 2 = 
B(B° -> 7T + 7r-)| tre e in units of 1(T 6 . 
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Figure 2: Fit to / dq 2 -j^(B° — > it l + ui) values obtained for three q 2 bins in Ref. 



In Fig. |21 we show our best fit to the CLEO data [5] for the B° — > tt~1 + V[ spectrum 
(in three q 2 bins). The data favor a rather lower value of F(0) than in our previous 
discussion jl], accounting for the lower magnitude of the tree amplitude in the present 
treatment. In Fig. El we show the comparison of the lattice data points with our best- 
fit form factor F + (q 2 ). As a consequence of the internal variations within the lattice 
results, a \ 2 °f about 5.5 (contributed by the lattice data) should be common for all fits; 
see Table II. Therefore, since the B — > ttIv spectrum is the second largest x 2 source, a 
significantly better overall fit can be achieved only if the measured B — > nlu branching 
ratios in the three q 2 bins are fitted better. This will require the addition of a fourth 
parameter to affect the shape of dB(B — > nlu)/dq 2 so that it is suppressed at both low 
and high q 2 ends and enhanced in the middle while relatively unchanged in the region 
13.6 GeV 2 < q 2 < 23.4 GeV 2 where lattice data exist; see Fig. |21 Consequently smaller 
tree amplitudes are implied and we regard them as disfavored by the lower limit as 
obtained earlier. 

V. HOW KNOWING THE TREE AMPLITUDE HELPS 

The ratio R n7T of the observed B° — > 7r + 7r~ branching ratio to its value in the presence 



7 



10 



15 



q 2 (GeV 2 ) 



20 



25 



Figure 3: Comparison of lattice data points with our best-fit form factor F + (q 2 ). Lat- 
tice data are from UKQCD (squares), APE (stars), Fermilab (circles) and JLQCD (di- 
amonds) . 



of the tree amplitude alone helps to establish the relative magnitude and strong and 
weak phase of the penguin amplitude in this process 0. On the basis of the previous 
determination of the tree amplitude [I] and the present world average for B(B° — > 7r + 7r~) 
we quoted [2E] Rmr = 0.62 ± 0.28, which indicated that tree-penguin interference was 
not required but, if present in the rate, would be destructive. The new information on 
\T\ allows us to refine this estimate to obtain = 0.87^o;2g, a value still consistent 
with both possibilities. 

The ratio \P/T\ of penguin to tree amplitudes quoted in Ref. j2£] was \P/T\ = 
0.28 ± 0.06. This ratio is useful in interpreting CP-violating asymmetries in the decay 
B° — > 7r + 7r~ (see, e.g., jHj). With the new world average [2E] B(B + — > K°tt + ) = 
(19.6 ± 1.4) x 10" 6 and the prescription \P/P'\ = (U/f K )X/(l - A 2 /2) we find for 
U = 130.7 MeV, f K = 159.8 MeV, and |2Zj A = 0.224 the values \P'\ = 4.28±0.16, \P\ = 
0.80 ±0.03, and \P/T\ = 0.35±g;g|. (Here the prime denotes a \AS\ = 1 amplitude.) The 
"penguin pollution" thus is slightly greater than estimated previously. Corrections to the 
CKM phase a obtained from the asymmetry parameter S n7T and the direct asymmetry 
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parameter A n7T both can be slightly larger than in Refs. 9\ and [26J. 

The tree/penguin ratio in B° — > K + n~ is also affected. By a similar analysis we 
found r = |T'/P'| = 0.173 ± 0.039 in Ref. |2S|; the new value is 0.151±g;g§f. A bound 
on the CKM phase 7 quoted in Ref. [2SJ relied on the lower limit of r, which is slightly 
raised, so the bound is strengthened slightly. Since it was only at the la level, we do 
not present it here. 

A further implication of the improved upper bound on T is a lower bound on C . 
Given the la bound |T + C\ > 2.89 based on the B + — > 7r + 7r° branching ratio (see Sec. 
Ill) and the la upper bound |T| < 2.65 based on the present analysis, we conclude that 
if C and T have a small relative phase ^U], then Re(C/T) ~ 0.1. 

VI. CONCLUSIONS 

The measurement of the B — > -niv spectrum by the CLEO Collaboration [S] has 
provided valuable information allowing us to improve the determination of the "tree" 
contribution to B° — > tt + tt~ . Combining this information with inclusive determinations 
of the CKM matrix element \V u b\ and lattice gauge theory calculations of the B —>■ n form 
factor F + (q 2 ), we have found B(B° — > 7r + 7r~)| troe = (5.25±o!o) x 10~ 6 , not significantly 
greater than the experimental value B(B° — > 7r + 7r~) = (4.55 ± 0.44) x 10~ 6 . The fit 
implies \V U b\ = (3.62 ± 0.34) x 10~ 3 , with an additional theoretical error of 15%. The 
relative strength of the penguin amplitude in this process, gauged using flavor SU(3) 
from the rate for B + — > K°tt + , is slightly larger than estimated previously, amounting 
to (35±§)% in amplitude. However, the need for strong destructive interference between 
this amplitude and the tree contribution is somewhat diminished in comparison with 
earlier estimates. 
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